function [DMRS_obs_mat,DMRS_Channel_H] = Function_Get_PSBCH_DMRS(NSLsymb, NFFT, cpLen0, cpLenR, NSLRB, l_DMRS, NSLID, input_signal)
 
	Msc_PSBCH = 72;                     % bandwidth of PSBCH in number of subcarriers
    NRBsc = 12;                         % Resource block size in the frequency domain, expressed as a number of subcarriers
    OFDM_Demod_Offset = 2;
    
    N_PRB = Msc_PSBCH / NRBsc;
    
	Hopping = 'Off';            % Hopping method ('Off', 'Group', 'Sequence')
    NCS = 0;
    
    % PSBCH subcarrier resources: central 72 subcarriers are allocated
    subixs_PSBCH = (0:1:Msc_PSBCH-1)' - Msc_PSBCH/2 +  NRBsc*NSLRB/2;
    
	Mode = 'psbch_V2X';
	[PSBCH_DMRS_Seq, DRSSymInfo] = Function_Generate_DMRS_Seq(N_PRB,  NSLID, NCS, Mode, Hopping);

    Output_signal = Function_Phy_OFDM_Demodulate_per_Subframe(NSLsymb, NFFT, cpLen0, cpLenR, NRBsc, NSLRB, input_signal, OFDM_Demod_Offset);
    
    DMRS_Seq_Rx = Function_Phy_Resources_Demapper(l_DMRS, subixs_PSBCH, Output_signal); % get the sequence (e.g. for broadcast: 144x1 [2 symbols of 72 subs vectorized)

    DMRS_obs_mat = -ones(length(DMRS_Seq_Rx)/length(l_DMRS),length(l_DMRS));
    DMRS_ideal_mat = -ones(size(DMRS_obs_mat));
    for lix = 1:size(DMRS_obs_mat,2)
        DMRS_obs_mat(:,lix) = DMRS_Seq_Rx((lix-1)*size(DMRS_obs_mat,1)+1:lix*size(DMRS_obs_mat,1));
        DMRS_ideal_mat(:,lix) = PSBCH_DMRS_Seq((lix-1)*size(DMRS_ideal_mat,1)+1:lix*size(DMRS_ideal_mat,1));
    end
    
    Channel_Est_Method = 'LS';
    [DMRS_Channel_H] = Function_Channel_Estimate_Process(DMRS_obs_mat,DMRS_ideal_mat,Channel_Est_Method);
    
end

    
    